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ABSTRACT 



Context. Red supergiant stars (RSGs) and yellow hypergiant stars (YHGs) are believed to be the high-mass counterparts of stars in the asymptotic 
giant branch (AGB) and early post-AGB phases. As such, they are scarcer and the properties and evolution of their envelopes are still poorly 
understood. 

Aims. We study the mass-loss in the post main-sequence evolution of massive stars, through the properties of their envelopes in the intermediate 
l—l, and warm gas layers. These are the regions where the acceleration of the gas takes place and the most recent mass-loss episodes can be seen. 
<J^ Methods. We used the HIFI instrument on-board the Herschel Space Observatory to observe sub-millimetre and far-infrared (FIR) transitions of 
p-h^ ■ CO, water, and their isotopologues in a sample of two RSGs (NMLCyg and Betelgeuse) and two YHGs (IRC+ 10420 and AFGL2343) stars. We 
>»✓ ' present an inventory of the detected lines and analyse the information revealed by their spectral profiles. A comparison of the line intensity and 
pH ] shape in various transitions is used to qualitatively derive a picture of the envelope physical structure. On the basis of the results presented in an 
f») ^ earlier study, we model the CO and 13 CO emission in IRC+ 10420 and compare it to a set of lines ranging from the millimetre to the FIR. 
I ■ Results. Red supergiants have stronger high-excitation lines than the YHGs, indicating that they harbour dense and hot inner shells contributing to 
O these transitions. Consequently, these high-./ lines in RSGs originate from acceleration layers that have not yet reached the circumstellar terminal 
j—j ] velocity and have narrower profiles than their flat-topped lower-./ counterparts. The YHGs tend to lack this inner component, in line with the picture 
C/2 , of detached, hollow envelopes derived from studies at longer wavelengths. NH3 is only detected in two sources (NML Cyg and IRC+ 10420), which 
• are also observed to be the strongest water-line emitters of the studied sample. In contrast, OH is detected in all sources and does not seem to 
' correlate with the water line intensities. We show that the IRC+ 10420 model derived solely from millimetre low-/ CO transitions is capable of 
reproducing the high-/ transitions when the temperature in the inner shell is simply lowered by about 30%. 

' Key words. Stars: AGB and post-AGB - supergiants - circumstellar matter - Submillimeter: stars 
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£^ . 1. Introduction that can be very dense. Mass-loss plays an important role in 

the evolution of these objects (de Jager 1998, Meynet & Maeder 
OO Red supergiant stars (RSGs) and yellow hypergiant stars (YHGs) 2003) . almost one half of their total initial mass can be ejected 

O . are thought to be stages in the post main-sequence evolution of during these late phaseSj affec ting in particular their possible 

stars with initial masses between -10 solar masses and 50 M later evolution into supernova events. Some objects have very 

*7J ; (e.g. de Jager 1998, Meynet & Maeder 2003, Levesque 2010). high mass _i oss rates , with episodic rates as high as M ~ 10" 3 

> . As such, RSGs and YHGs are the massive counterparts of AGB M& ^ wh0e the circumst ellar envelopes around others are 

and (early) post-AGB stars. However, the very different prop- very diffus6i coll - eS p nding to rates under 10" 7 M yr _1 (e.g. 

erties of evolved massive stellar objects (RSGs, YHGs, Wolf- C astro-Carrizo et al. 2007 - CC07 thereafter, Quintana-Lacaci 

Rayet stars, luminous blue variable stars, Cepheid-hke variables, 2008i Mauron & Josselin 201 1). From the theoretical point of 

supernovae etc), particularly their distribution in the H-R dia- vieWi yery high rates ^ also exp ected, at least episodically, ow- 

gram, is hard to interpret with a simple description of the evolu- ing t0 a comb ination of high radiation pressure and atmospheric 

tion in this phase. activity (e.g. Josselin & Plez 2007; as happens in the AGB) or 

Both RSGs and YHGs are known to show complex mass- t o the intrinsic instabilities characteristic of the yellow phases 

loss phenomena, which form circumstellar envelopes (CSEs) G f t h e late evolution of high-mass stars (the so-called "yellow 

void"; see e.g. de Jager 1998). 
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* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- Molecular lines have been one of the most powerful tools for 
tant participation from NASA. the study of the properties of CSEs around AGB and post-AGB 
** Figures 10-13 and Appendix are only available in electronic form stars. Similar data on RSGs and YHGs are rare, but in some 
at http://www.edpsciences.org cases molecular lines have yielded a quite comprehensive insight 
*** The spectra presented in this paper are available in electronic form into their CSEs. For instance, using model fitting of mm- and 
at the CDS submm-observations of several CO rotational lines in the RSG 
1 Deceased 14 January 201 1 VY CMa, Decin et al. (2006) deduced mass-loss rates as high as 
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3xl0~ 4 M© yr -1 , that varied on timescales of about 1000 yr. VY 
CMa is known to harbour a wide variety of molecular species, 
illustrating the very complex chemistry at play in its CSE (e.g. 
Tenenbaum et al. 2010 and references therein). Other less well- 
studied RSGs, such as VX Sgr and NML Cyg, also show high 
mass-loss rates of the order of ~ 10~ 4 M© yr~'(e.g. De Becket al. 
2010). An example of a low-M RSG is Betelgeuse (aOri), with 
a value of about 10 7 -10 6 M Q yr -1 , although this value is uncer- 
tain because its molecular emission may be weaker than in other 
objects owing to its molecular underabundance (see Huggins et 
al. 1994, Mauron & Josselin 201 1). The YHGs IRC+10420 and 
AFGL 2343 have also been well-studied in molecular emission. 
CC07 (see also Quintana-Lacaci et al. 2008) compared line- 
emission models with high-resolution maps of the CO 7=1-0 
and J -2-1 lines, deriving large M variations on timescales of 
about 1000 yr and M maxima in excess of 10~ 3 M Q yr -1 . These 
objects also have many molecular lines and a very rich chemistry 
(Quintana-Lacaci et al. 2007). The total masses in the molecule- 
rich shells are particularly high for these two YHGs, M tot i 1 
M Q . The total masses of the molecular envelopes around RSGs 
are somewhat lower, between 0. 1 solar masses and 1 M© for VY 
CMa, VX Sgr, and NML Cyg. 

Despite the observational progress that has been made to 
date, the available data still fail to provide information on some 
basic parameters. Low-/ CO lines are good tracers of both the 
circumstellar mass distribution and kinematics, but these easily 
excited lines cannot probe warm gas with temperatures > 100 
K, which are expected to be present and sometimes dominant in 
these shells. The temperature itself is not well-determined un- 
der these conditions, which may affect the determination of the 
mass-loss rate and the total mass. To properly study these warm 
components, it is therefore necessary to observe lines in the far- 
infrared (FIR), involving level excitations comparable (in tem- 
perature units) to these moderate kinetic temperatures. 

This paper presents Herschel/HIFI observations in the FIR of 
molecular lines from two RSGs, Betelgeuse and NML Cyg, and 
two YHGs, IRC+10420 and AFGL 2343. The observations are 
part of the Herschel guaranteed time key program HIFISTARS 
(PI V. Bujarrabal), devoted to the study of high-excitation molec- 
ular lines in (low- and high-mass) evolved stars. Companion 
observations of VY CMa, which are particularly rich and com- 
plex, will be discussed in another paper (Alcolea et al., in prep.). 
We present data of the 7=6-5, 7=10-9, and 7=16-15 lines of 
12 CO and 13 CO, which are good probes of the excitation. We 
also discuss our observations of other molecules, including wa- 
ter vapour, OH, and NH3, which are particularly useful to under- 
standing more clearly the chemistry in these unusual objects. 

2. Observations and data processing 

The observations presented in this paper were obtained with the 
Heterodyne Instrument for the Far Infrared (thereafter HIFI, de 
Graauw et al. 2010), on-board the Herschel Space Observatory 
(Pilbratt et al. 201 1). Similarly to other observations presented in 
the framework of the HIFISTARS key program (e.g. Bujarrabal 
et al. 2012), the main target lines of this project were collected 
through a handful of HIFI settings. Thanks to the instantaneous 
spectral coverage of 2.4 - 4 GHz offered by the Wide-Band 
Spectrometer (WBS), and since HIFI works with double-side- 
band (DSB) mixers, several other bonus lines could be covered 
simultaneously. A summary of the settings used and their main 
characteristics is given in TableQ] 

All data were taken in the double-beam switching (DBS) 
mode, consisting of chopping and nodding between the source 



Table 1. Summary of observational settings and telescope char- 
acteristics. The setting names are provided for cross-reference 
with other HIFISTARS results, and to ease the association of 
data with a given Herschel observation identifier (obsid). These 
obsid's are given in the figures of the appendix. For each set- 
ting, we indicate the corresponding HIFI band, the local oscil- 
lator frequency tuned for the observation and the system noise 
temperature achieved on average, together with the typical beam 
properties applying to each of them. 



Setting 
name 


LO freq. 
(GHz) 


Typical T sys 
(K, DSB) 


Band 


Beam 
size 


5 e ff 


14g 


564.56 


93 


IB 


37.5" 


0.754 


12f 


653.55 


131 


2A 


32.4" 


0.752 


17 


686.42 


142 


2A 


30.9" 


0.745 


07e 


1106.90 


416 


4B 


19.1" 


0.742 


06d 


1157.67 


900 


5A 


18.3" 


0.639 


05c 


1200.90 


1015 


5A 


17.6" 


0.633 


03b 


1757.68 


1580 


7A 


12.1" 


0.740 


16 


1838.31 


1300 


7B 


11.5" 


0.736 



line-of-sight and blank sky positions 3 arcmin on either sides of 
the source. This mode allows very fast switches between the sky 
positions, and therefore mitigates the detector drifts. Residual 
baseline ripples can however remain, especially in the HIFI 
bands 6 and 7 where the detectors used (hot electron bolome- 
ters, hereafter HEB) are more susceptible to instabilities. 

We used the standard pipeline products that had been gener- 
ated and made available at the Herschel Science Archive (HSA). 
These could be accessed via the HIPE softwarfl The HIFI 
pipeline brings all data onto a single-side-band 7^ intensity 
scale. A side-band ratio of unity was assumed for all the lines 
reported in this paper. All individual spectra from the so-called 
levell were checked for outliers, and averaged for each of the 
two orthogonal mixer polarisations available on HIFI. In some 
rare cases, one of the polarisations was useless owing to an un- 
pumped mixer, so its data had to be discarded. The sky positions 
observed by the two polarization mixers are not strictly the same, 
the misalignment varying between ~ 6" in the lower bands to 
1" in the higher bands (Roelfsema et al. 2012). Any such mis- 
alignment could in practice manifest itself (among other cali- 
bration effects) as intensity imbalance between the two polarisa- 
tions. We however averaged the two spectra whenever possible 
to improve the sensitivity. In practice, this implies that the ef- 
fective beam size of the averaged data is slightly larger than the 
nominal spatial resolution. 

A special treatment had to be applied to some of the data 
taken in the HEB bands owing to the standing wave affecting 
in particular these detectors. This spectral artefact is not opti- 
cal in nature, hence cannot be optimally treated by applying the 
standard defringing methods to sine-wave baseline distortion. A 
novel method is currently being developed to perform an alterna- 
tive pipeline correction of these features (Higgins 201 1), which 
may, however, be ineffective when the line widths become simi- 
lar to those of the artefact structure. Because this particular situ- 
ation applies to the sources treated in this study, we decided not 
to adopt this approach. We developed instead a semi-automated 
procedure that discards individual spectra where the contribution 
of the baseline ripple is significantly stronger than the expected 
radiometric noise. This results in the rejection of up to half of 



1 HIPE is a joint development by the Herschel Science Ground 
Segment Consortium, consisting of ESA, the NASA Herschel Science 
Center, and the HIFI, PACS, and SPIRE consortia. 
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T mb (K) NML Cyg - 12 C0 and HCN T mb (K) NML Cyg - 13 C0 and NH 3 




-50 50 -50 50 



Velocity (km/s) Velocity (km/s) 

Fig. 1. Individual line spectra collected with Herschel/HIFI on 
NML Cyg for CO, 13 CO, NH 3 , and HCN. The dashed line indi- 
cates the source systemic velocity. 



the non-averaged spectra, and therefore increases the noise in 
the final data, in spite of the more reliable baseline. 

Finally, all spectra were converted onto a r m b scale, by 
applying the main beam efficiencies reported by the HIFI 
Instrument Control Centre, which are reproduced in TableQ] We 
based our absolute-calibration accuracy estimates on the error 
budget reported by Roelfsema et al. (2012). On top of this, we 
considered an additional contribution arising from the standing 
wave described above for bands 6 and 7. Adding up all eiTor con- 
tributions in quadrature only provided lower limits to the typical 
calibration errors. In practice, the calibration error that we as- 
sumed in this study is 15% in bands 1 and 2, 20% in bands 3 
to 5, and 30% in bands 6 and 7. The collection of all full-band 
spectra is compiled in Appendix lAl 

3. Results 

Table [2] summarises the main stellar parameters of the sample 
studied here. In the following sub-sections, we give an overall 
inventory of all probed and detected species, together with the 
extracted line characteristics in the form of the line peak and 
integrated intensities (see Tables [3] and [4] for a compilation of 
these results). The high spectral resolution of HIFI also allows 
us to distinguish and discuss the contribution of particular shell 
layers to the emission of the respective species. 

3.1. NML Cyg 

NML Cyg is a red supergiant of spectral type M6. It has an in- 
termediate mass-loss rate of order 8.7xlO~ 5 M yr~'(De Beck et 
al. 2010). 

NML Cyg exhibits the richest spectrum of all four sources 
presented here, although, among the full sample of RSGs, it is 
largely surpassed by the line density of VY CMa (Alcolea et al., 
in prep.). Apart from this latter, it is also the source with the 



T mb (K) NML Cyg - o-H a and OH T mb (K) NML Cyg - p-H s 




-50 50 -50 50 



Velocity (km/s) Velocity (km/s) 

Fig. 2. Same as Fig. Q] for OH, H2O and isotopologues in 
NML Cyg 

T mb (K) NML Cyg - SiO and SO T mb (K) NML Cyg - SiO and SO 




-50 50 -50 50 



Velocity (km/s) Velocity (km/s) 

Fig. 3. Same as Fig. [2] for SiO, 29 SiO, 30 SiO and SO in 
NML Cyg. 
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Table 2. Observed coordinates and additional stellar parameters for the two RSGs and the two YHGs considered in this study: 
spectral type, distance and luminosity. References are: 1: Mauron & Josselin (2011); 2: Schuster et al. (2009); 3: de Jager (1998); 
4: Hawkins et al. (1995), Reddy & Hrivnak (1999). 



Source a (J2000) 


d (J2000) 


Spectral type 


Distance (pc) 


Luminosity (L G ) 


References 


Red supergiants 


Betelgeuse 05 55 10.3 


+07 24 25.4 


M2 


150 


4 10 i 


1 


NMLCyg 20 46 25.5 


+40 06 59.6 


M6 


1700 


6 10 5 


2 


Yellow hypergiants 


IRC+ 10420 19 26 48.0 


+11 21 16.7 


F8 


5000 


7 10 5 


3 


AFGL2343 19 13 58.6 


+00 07 31.9 


G5 


6000 


6 10 5 


4 



strongest line emission when considering species and transitions 
individually (see Tables[3]and|4] and Fig.[9]l. While the detection 
of water vapour lines in other sources is mostly limited to the 
ground-state transitions of both ortho- and para-^O, NMLCyg 
shows transitions with energy levels in excess of 2000 K. One 
noticeable detection is that of the H2O maser line at 658 GHz. 
The measured flux of ~400Jy agrees well with the first detec- 
tion of this line by Menten & Young (1995), indicating that 
there is no evidence of large variability. Water maser emission 
at 22 GHz has also been reported in this source (e.g. Nagayama 
et al. 2008). While the velocity distribution of this latter exhibits 
two shifted components at -25 km s and +5 km s respec- 
tively, the 658 GHz maser is peaked at the star velocity of Vlsr = 
km s _1 , with a relatively narrow profile compared to the other 
submm transitions, and a blue-shifted shoulder. This is in con- 
trast to the rest of the water vapour lines, where the blue side 
is usually narrower than the red side. The narrower line width 
shows that the maser originates from inner layers still in accel- 
eration. 

NML Cyg also displays a uniquely rich spectrum of SO, ro- 
tational transitions from the ground, and vibrationally excited 
states of SiO, as well as isotopes of the SiO, though only from 
the ground vibrational state (see Fig. [3}. The detection of high-/ 
transitions of SiO confirms the presence of warm and dense gas 
at the centre of the envelope. This is a noticeable observation 
since, as shown in the following sections, the other sources seem 
to have rather hollow envelopes. This result is corroborated by 
the presence of very high excitation H2O transitions, and tells us 
that the mass-loss rate is probably still very high today. Finally, 
NML Cyg is one of the few sources showing a strong emission 
in NH 3 , and the only RSG together with VY CMa to exhibit de- 
tectable emission from other hydrides such as HCN and H2S. 

The velocity structure of the various species and transitions 
displays a great variety of profile shapes. As is well-known, the 
line shape essentially depends on the line opacity, but also on the 
velocity profile of the layers where the line is emitting. Typically, 
the transition from a rectangular (flat-top) to a parabolic profile 
reflects the effect of optical depth in a constant velocity flow, 
while the transition to a triangular, narrower profile comes about 
when one probes the inner and warmer acceleration region with 
higher excitation lines. A recent study of submm molecular lines 
of CO and H 2 in O-rich AGB stars (Justtanont et al. 2012) in- 
dicates that the line width tends to be smaller for high excitation 
lines formed closer to the star, pointing to an accelerated outflow. 

In NMLCyg, this kind of transition is particularly pro- 
nounced when looking at the evolution between the 7=6-5, 
7=10-9, and 7=16-15 transitions of 12 CO and 13 CO, respec- 
tively (see Fig. [TJ. That the same trend is observed in both op- 
tically thick and thin isotopes suggests that this is not primarily 



an opacity effect and we indeed detect the velocity gradient in 
the innermost layers of the envelope. The same behaviour is ob- 
served for increasingly higher energy levels of the /7-H2O and 
p-H^O transitions (see Fig. |2|. Opacity, however, also plays 
a role in the line shape, as can be seen from the narrowing 
of the line profile between the respective 7=10-9 transitions 
of 12 CO and 13 CO, which have similar excitation temperatures. 
Interestingly, the 1834 GHz OH line profile closely resembles 
that of the mid-7 CO transitions, suggesting that its emission 
arises predominantly in the outer layers where the terminal ve- 
locity has been reached. 

Another peculiarity of the velocity structure can be observed 
in the NH3 line profile (see Fig.[TJ). Unlike the other lines studied, 
that exhibit an overall single-peaked profile, the NH3 emission 
has two separate peaks on either side of the systemic velocity. In 
particular, this structure provides evidence for a component ly- 
ing at velocity around -18 km s . We note that this blue-shifted 
velocity component coincides with a prominent contribution of 
the v=l, 7=2-1, SiO maser emission at 43 GHz (Boboltz & 
Marvel 2000), showing up at sub-arcsec scales as a marked bi- 
polar structure also reported in mid-IR maps of the dust emission 
(e.g. Monnier et al. 2004). Menten et al. (2010) already empha- 
sized the peculiar behaviour of ammonia in other giant stars such 
as VYCMa or IRC+ 10420, reporting unexpectedly high abun- 
dances in their envelope. We speculate that the bulk of the NH3 
abundance could actually be quite widely distributed over the 
circumstellar shell. 

Finally, we note the various ISM absorption signatures seen 
only in the two ground-state transitions of 0-H2O and /7-H2O. 
While these features were also observed in low-7 CO transitions 
(Kemper et al., 2003), they are not present in the other transitions 
probed by Herschel, confirming that they arise from either a dif- 
fuse or cold gas component. We also note that there is a narrow 
emission component in the ground state 0-H2O and P-H2O at 
5 km s _1 . This small peak was also observed in the 12 CO 7=2-1 
line of Kemper et al. (2003) and probably originates from a small 
ISM clump along the line-of-sight. This interpretation is sup- 
ported by recent distance measurements of NML Cyg (Zhang et 
al. 2012) and foreground interstellar gas (Rygl et al. 2012). 

3.2. Betelgeuse 

Betelgeuse (a Ori) is a red supergiant of spectral type M2. It has 
a relatively low mass-loss rate of order 2x 1O~ 7 M yr _1 (De Beck 
et al. 2010). 

Betelgeuse appears as a source that is particularly strong in 
I2 CO and 13 CO emission. In absolute intensities, it is the bright- 
est of all targets reported here for these molecules. In contrast, 
it exhibits faint H2O and OH lines, especially in relation to the 
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Table 3. Spectral line results for NMLCyg and IRC+ 10420. For each line we indicate the corresponding upper energy level and 
rest frequency, together with the measured line peak intensity, line integrated intensity, and the velocity range used to compute this 
latter. 



Species 


Transition 


£ f 

^upp 

(K) 


Rest freq. 
(GHz) 


Peak <7) 
(mK) 


NML Cyg 
Integ. intensity 
(Kkm s -1 ) 


Vel. range <2) 
LSR (km s _1 ) 


Peak <7) 
(mK) 


IRC+ 10420 
Integ. intensity 
(Kkm s~') 


Vel. range' 2 ' 
LSR (km s~') 


u CO 


7=6-5 


116 


691.473 


471(iJ) w 


19.6 


[-34; 33] 


364(20) (3) 


19.1 


[37; 115] 




7=10-9 


304 


1151.985 


479(^9) <4) 


19.1 


[-33; 30] 


243(44) (4) 


14.7 


[36; 114] 




7=16-15 


752 


1841.346 


666(47) <4) 


17.6 


[-38; 28] 


133(44) (5) 


6.57 


[38; 111] 


"CO 


7=6-5 


1 1 1 


661.067 


95(7) (4) 


4.06 


[-34; 33] 


57(6) (4) 


2.79 


[36; 113] 




7=10-9 


291 


1101.350 


141(i7) (4) 


3.80 


[-25; 31] 


60(20) <4) 


2.8 


[43; 107] 




7=16-15 


719 


1760.486 


135(<«) (6) 


4.75 


[-31; 27] 


<138<'- 5 > 


- 


- 


o-H 2 


7jr a ,jr c = l i,o-l 0,1 


27 


556.936 


386(6) <3) 


12.2 


[-23; 37] 


168(6) <3) 


6.80 


[43; 118] 




J K, K c —3, 2 — 22 1 


215 


1153.127 


848(4i) <4) 


28.8 


[-23; 35] 


333(42)< 4) 


13.7 


[38; 111] 




7k„,a: 1 ; = 32,i-3i,2 


271 


1162.912 


869(42) <4) 


29.7 


[-31; 34] 


195(4i) (4) 


8.41 


[40 ;112] 


p-H 2 


7jr a ,*i=li,i-Oo,o 


53 


1113.343 


993(20 (3) 


30.8 


[-23; 38] 


329(2i) (3) 


12.7 


[46; 115] 




J k" v — 4t t — t-1 -\ 


454 


1207.639 


726(47) <5) 


22.9 


[-32; 34] 


<144 (1,4) 










952 


1762.043 


526(45) (4) 


14.0 


[-29; 30] 


<180 (1 - 4) 


- 


- 


O-H2O V2=l 


J k k~ — 1 1 n~ln 1 


2326 


658.004 


882(9) <3) 


9.53 


[-26; 18] 


<15 (1 - 5) 






p-H 2 V2=l 


7y a ,x e =li,i-Oo,o 


2352 


1205.788 


112(44) <4) 


2.64 


[-18; 17] 


< 144<i,4) 






p-H' 8 


7x a ,x c =li,i-Oo,o 


53 


1101.698 


185(iS) (4) 


5.05 


[-24; 37] 


<60 <14) 


- 


- 


SiO 


7=16-15 


283 


694.294 


118(i4) (4) 


3.42 


[-32; 30] 


<35 (1 ' 5 » 


- 


- 


SiO v=l 


7=13-12 


1957 


560.325 


17(3) (S) 


0.39 


[-25; 27] 


<15 u>4 ' 




— 




7=15-14 


2017 


646.429 


30(5) (5) 


0.50 


[-25; 25] 


1 tr(\ SI 

<15 u,: " 






29 SiO 


7=13-12 


187 


557.179 


30(4) (5) 


0.95 


[-25; 26] 


<15<'> 4 > 




- 




7=26-25 


722 


1112.798 


50(i7) (5) 


1.10 


[-21; 13] 


<60 <14) 






30 SiO 


7=26-25 


713 


1099.708 


59(i4) (5) 


0.77 


[-17; 24] 


<60 <14) 






SO 


J K = 13 14 — I2j3 


193 


560.178 


13(4) (5) 


0.40 


[-22; 23] 


<15 (1 - 4) 








7^=13,2-12,; 


194 


558.087 


13(4) (5) 


0.56 


[-35; 29] 


<15<»- 4 ) 


- 


- 




7 K =13, 3 -12,2 


201 


559.320 


20(5) (5) 


0.49 


[-22; 25] 


<15<'> 4 > 








7 K =15, 6 -14,5 


253 


645.875 


21(4) (5) 


0.54 


[-21; 22] 


<15<»- 5 > 


- 


- 


H 2 S 


7^,^=31,2-22,1 


137 


1196.012 


149(33) m 


3.83 


[-13; 22] 


<144 <1.4) 






HCN 


7=13-12 


386 


1151.452 


86(27) (6) 


3.48 


[-35; 32] 


<126 <L4) 






H 13 CN 


7=8-7 


147 


690.552 


<18 (1,6) 






<35 <1,5) 






NH 3 


7k=1o-Oo 


27 


572.498 


102(3) (4) 


4.08 


[-38; 38] 


93(5) <4) 


2.93 


[41; 116] 


OH 2 n 1/2 


7=3/2-1/2 


270 


1834.747 


800(46) <4) 


31.4 


[-25; 36] 


599(5S) <4) 


33.5 


[37; 113] 



For o-H 2 0, the energies are relative to that of the lowest ortho level (7^,^ = 10,1), since 0-H2O and /7-H2O behave as different species. 
' Non-detections are given as 3-cr, <2) Range over which the integrated intensity is computed, <3) Computed in bins of order 0.5 km s~' , 
' Computed in bins of order 1 km s~' , <5) Computed in bins of order 2 km s~' , <6) Computed in bins of order 3 km s~' , 
' Computed as the line maximum minus 1.5x the noise rms, indicated in italics between brackets (l-cr) 



CO line intensity (Fig.|4j see also Fig. [9}. Betelgeuse has a par- 
ticularly small distance, 150 pc, which, despite its low mass- 
loss rate, explains its high CO line intensities. Owing to this 
weak mass loss, the envelope density is lower, and the weak 
water emission could be a consequence of enhanced photo- 
dissociation by stellar and interstellar ultra-violet (UV) photons 
as H2O is more easily photodissociated than CO (e.g. Justtanont 
et al. 1999). We recall that this source has very strong fine- 



structure lines of both neutral atoms and low-energy ions (e.g. 
Haas & Glassgold 1993, Castro-Carrizo et al. 2001), which are 
thought to be abundant in the inner envelope owing to the UV 
emission of the extended chromosphere of the star. 

As in NML Cyg, the 12 CO and 13 CO line profiles become 
narrower as the excitation temperature increases, which again 
suggests that the high-transition lines form predominantly in the 
acceleration region. Betelgeuse is thought to have a particularly 
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Table 4. Same as Table[3]for the spectral line observed in AFGL2343 and Betelgeuse 



Betelgeuse AFGL2343 



Species 


Transition 




E' 

^upp 

(K) 


Rest freq. 
(GHz) 


Peak (6) 
(mK) 


Integ. intensity 
(Kkms- 1 ) 


Vel. range*^ 
LSR (km s" 1 ) 


Peak (6) 
(mK) 


Integ. intensity 

(Kkms- 1 ) 


Vel. range^' 
LSR (kms 1 ) 


12 CO 


.7=6-5 




116 


691.473 


654(7 7)^ 


13.2 


[-12; 24] 


261(iJ) w 


11.9 


[61; 135] 




j = 10-9 




304 


1151 985 


762(54) (3) 


13.4 


[-12- 211 


21 1C2iTt (5) 


6.94 


[61; 137] 




7=16-15 




752 


1841.346 


866(5S) (3) 


13.6 


[-11; 22] 


< 111< 15 > 






13 co 


7=6-5 




111 


661.067 


136(7) (3) 


2.83 


[-10; 18] 


87(6) (4) 


3.00 


[62; 133] 




7=10-9 




291 


1 1 A 1 OCA 

1 101.350 


iov(2uy ' 


2.45 


[-9; 23] 


< 156 u ' ' 








7=16-15 




719 


1760.486 


187(50) (5) 


2.68 


[-9; 22] 


< 105 (1 - 5) 






o-H 2 




lo,i 


27 


556.936 


35(4)< 4 > 


0.61 


[-6; 21] 


36(4)< 4 > 


1.41 


[62; 132] 




Jk^.Kc =32,1- 


3l.2 


271 


1162.912 


< IV2PS> 






98(iJ) (5) 


4.26 


[73 ;131] 


p-H 2 


■^AT a ,*:c = 1 l.l" 


Oo.o 


53 


1113.343 


72(i6~) (4) 


1.54 


[-9; 21] 


49(i3) (4) 


1.93 


[77; 130] 


OH 2 n 1/2 


7=3/2-1/2 




270 


1834.747 


208(i4) (4) 


4.45 


[-10; 20] 


362(i0) (5) 


14.0 


[64; 129] 



For o-HtO, the energies are relative to that of the lowest ortho level (Jk i1 jc c = Io,i), since 0-H2O and /7-H2O behave as different species. 
' Non-detections are given as 3-cr, <2) range over which the integrated intensity is computed, (3) Computed in bins of order 0.5 km s -1 , 

Computed in bins of order 1 km s , <5) Computed in bins of order 2 km s , 
'* Computed as the line maximum minus 1.5x the noise rms, indicated in italics between brackets (1-cr) 




Fig 



Velocity (km/s) 

4. Same as Fig.|2]for lines detected in Betelgeuse. 



hot inner circumstellar region (Rodgers & Glassgold 1 99 1 ), with 
temperatures higher than about 1000 K out to 10 15 cm from the 
star (i.e. occupying about 1"). The evolution in the line profile 
can be clearly seen if one compares the spectra reported here to 
the nearly flat-top shape of the 12 CO 7=3-2 spectrum observed 
by Kemper et al. (2003). There are also some small spectral com- 
ponents in our CO profiles sitting on top of the overall parabolic 
or triangular profiles. Those small structures could be related 



to clumps in the inner envelope, which have been reported in 
mid-IR observations of the inner 3"of the nebula (Kervella et al. 
2011). 

3.3. IRC+10420 

IRC+10420 is a yellow hypergiant, of spectral type F8. It has 
a strongly variable mass-loss rate, with highest values around 
3xlO- 4 M yrHCCO?). 

IRC+10420 presents the second richest species inventory 
after NMLCyg (Fig. [5]), although it features none of the SO, 
SiO, or C-bearing molecules detected there (other than CO and 
13 CO). Noticeably, strong lines of NH3 and OH are observed, 
where hydroxil is the brightest of all detected lines in this source. 
Only a handful of water lines are detected on top of the funda- 
mental levels of 0-H2O and p-H 2 (see Table EJ. While the CO 
lines have roughly parabolic line profiles, the rest of the detected 
species display a strong asymmetry in the form of a blue-shifted 
intensity excess. This feature is particularly clear in the H2O 
and NH3 profiles (Fig. [5] see also Menten et al. 2010), but only 
mildly observed in OH. A similar trend was already reported 
in mm-wavelengths transitions of e.g. CS, H I3 CN, or 29 SiO by 
Quintana-Lacaci et al. (2007). CC07 demonstrated that the main 
contribution of this spectral component arises from a region in 
the western halo of the envelope (see their Fig. 8). This isolated 
shell component was also reported in high-resolution mid-IR 
maps by Lagadec et al. (201 1 ). We propose that this spectral fea- 
ture is due to an enhanced abundance of certain species in these 
layers, as the model presented by CC07 predicts only moderate 
(~50K) temperatures at this radius. H2O and NH3 are often be- 
lieved to appear as a result of shock-induced chemistry, suggest- 
ing that this gas layer could be composed of shocked material. 

Section|4]presents a revised modelling of the whole set of CO 
data and discusses the derived structure of the envelope around 
IRC+10420. 

3.4. AFGL 2343 

There is some controversy about the nature of AFGL 2343, 
which has been proposed by certain authors to be a post-AGB 
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T mb (K) IRC+10420 - H 2 0, NH a and OH T mb (K) IRC+10420 - 12 C0 and 13 C0 




50 100 50 100 



Velocity (km/s) Velocity (km/s) 

Fig.5. Same as Fig.|2]for lines detected in IRC+10420. 



T mb (K) AFGL 2343 - H a and OH T mb (K) AFGL 2343 - 12 C0 and 13 CO 




50 100 150 



Velocity (km/s) 

Fig. 6. Same as Fig.|2]for lines detected in AFGL 2343. 



(low-mass) object (see Ferguson & Ueta 2010 and references 
therein), located at a relatively short distance (~ 1 kpc, instead 
of ~ 6 kpc if it is a luminous hypergiant; note that no reli- 
able parallax measurement of this source exists). The majority 
of published papers favour the massive-star hypothesis, in par- 
ticular those reporting studies of the circumstellar shell around 
AFGL 2343 (Reddy & Hrivnak 1999, Gledhill et al. 2001, CC07, 
Quintana-Lacaci et al. 2008), since the shell seems very ex- 
tended and similar to that of the well-studied IRC+10420. The 
fast and mostly isotropical expansion of its envelope also sup- 
ports the YHG nature of AGL2343. Accordingly, we also as- 
sume that AFGL 2343 is a YHG, of spectral type G5. CC07 re- 
ported a variable and very high mass-loss rate, with values as 
high as 3xl0~ 3 M G yr _1 in the past, but having shown a sharp 
drop to a still considerable ~ 4xlO~ 5 M yr~' about 1200 years 
ago. 

AFGL 2343 exhibits the faintest set of submm lines of the 
sample studied here. In particular, it is the only source where no 
emission is detected in the 12 CO 7=16-15 and 13 CO 7=10-9 
transitions. This is in line with the modelling results reported 
by CC07, where AFGL 2343 appears as a detached shell, with 
relatively low temperatures in general and low density in the 
central layers. A similar picture is derived from high resolution 
mid-IR maps (Lagadec et al. 2011). CC07 deduced a temper- 
ature of lower than about 20 K for the densest shell and be- 
tween ~300 K and 30 K for the much-less-dense inner compo- 
nent. These low temperatures may however be underestimated 
given the detection of several transitions with energy levels near 
-300 K (e.g. 12 CO 7=10-9, o-H z O 3 2 ,i-3 u or OH). The lack of 
dense enough warm layers in this source is corroborated by the 
12 CO 7=10-9/7=6-5 line and integrated line ratios (-0.58 and 
~0.5, respectively). These values are significantly smaller than 
the corresponding line and integrated line ratios for IRC+10420, 
which reach ~0.77 and ~0.7 1 respectively, confirming that most 
of the circumstellar gas in AFGL 2343 is cooler than the domi- 
nant circumstellar shells in IRC+10420. 

The various line profiles have a red-shifted excess at Vlsr 
~ 125 km s _1 . The contrast between this narrow component and 
the emission around the star systemic velocity is higher for 
higher excitation lines. The same behaviour was observed in 
the mm-wave lines of several molecules (e.g. SiO, Quintana- 
Lacaci et al. 2007). In contrast, this excess tends to disappear in 
lower excitation lines and is for example hardly seen in the 12 CO 
7=1-0 and 7=2-1 spectra (Bujarrabalet al. 2001). Interestingly, 
the 12 CO 7=10-9/7=6-5 line ratio at this velocity increases to 
~0.81, implying that there is a gas component in AFGL 2343 
with excitation conditions at least comparable to those present 
in IRC+10420. 

The origin of this spectral component remains unclear. In the 
interferometric maps of 12 CO 7=1-0 and 7=2-1 by CC07, this 
feature appears as a compact clump located in the same direc- 
tion as the central star and moving almost at the extreme pos- 
itive LSR velocity. It would then correspond to a part of the 
envelope placed basically behind the star. Quintana-Lacaci et 
al. (2008) studied this component in more detail using maps of 
29 SiO 7=2-1 and HCN 7=1-0 emission. They concluded that 
this feature comes from a component at relatively small radius 
of high density and temperature, and deduced it to be a very 
clumpy, thin shell. These authors suggest that the high-excitation 
condensations probably result from anisotropic shock propaga- 
tion, particularly by analogy with the case of IRC+10420, in 
which signs of shocks running within the envelope were also 
found from SiO emission maps. 
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We speculate that the spectral feature found at 125 km s 
could also result from a shock between a hemisphere of the cir- 
cumstellar shell and nearby interstellar gas. It is remarkable that 
the LSR velocity of AFGL 2343 is positive and high in absolute 
value, ~100km s _1 , so the object is receding from us at high ve- 
locity. On the other hand, the low-/ CO observations (Bujarrabal 
et al. 2001) show that there is interstellar emission at moder- 
ate LSR velocities, between 5 km s _1 and 10km s _1 . We cannot 
be sure whether this gas is really close to the YHG, but if this 
were the case, the strong difference in the projected velocities 
should produce a shock in the shell hemisphere placed beyond 
the star, which would show an LSR velocity equal to approxi- 
mately the stellar velocity plus the circumstellar expansion ve- 
locity, i.e. about 125 km s _1 . 

4. Modelling of the CO emission in IRC+10420 

4.1. The model 

CC07 developed a detailed model of the structure and kinemat- 
ics of the CSE around IRC+10420 based on interferometric CO 
7=1-0 and / =2-1 maps, together with the sub-mm observa- 
tions of CO transitions up to 7=6-5 performed by Teyssier et 
al. (2006). 

We performed numerical calculations of CO line excitation 
and radiative transfer using a code very similar to that described 
in CC07. Despite the evidence of departure from spherical sym- 
metry discussed in CC07 and Section 13.31 this model assumed 
isotropic mass loss, and therefore used a one-dimensional code. 
The line excitation and the population of the rovibrational lev- 
els were calculated using an large-velocity-gradient (LVG) ap- 
proach, which is justified in view of the large macroscopic ve- 
locities found in this object. The brightness distribution for each 
line was calculated by solving the 'exact' radiative transfer equa- 
tions, so radiative interaction between distant points and opacity 
effects were fully taken into account. The brightness distribution 
was later convolved with the corresponding telescope beams, 
which were assumed to be Gaussian. 

The original model by CC07 consisted of two detached 
shells, each one corresponding to large mass loss episodes 
spaced at interval of 1200 years. It also assumed that the mass 
loss ceased 200 years ago. These respective mass ejections ex- 
pand at high velocities (37 km s _1 and 25 km s for the inner 
and outer shells, respectively), as expected given the high lu- 
minosity of these objects. Following CC07, the inner shell is 
assumed to be located between 2.5xl0 16 cm and 1.2xl0 17 cm 
and have a mass-loss rate of 3xl0~ 4 M Q yr _1 . The outer shell is 
located between 2.2xl0 17 cm and 5.2xl0 17 cm and has a rate 
of 1.2xl0~ 4 M Q yr _1 . The kinetic temperatures were described 
by a power law 7\(r) = rn(r/10 17 cm) Q ', where Tyi is the tem- 
perature at 10 17 cm from the star. The gas temperatures derived 
by CC07 were high, T l7 = 230 K and T Y1 = 100 K for the in- 
ner and outer shells, respectively. Values of a 1 were found. 

However, we note that low-/ transitions are not good tracers of 
high-excitation regions, so that HIFI observations are necessary 
to accurately constrain the temperature. 

4.2. Model results: Comparison with Herschel/HIFI data 

We compared the outcome of high-/ line profiles predicted by 
this CSE model with our Herschel/HIFI data. Although that 
model was derived from data of low-excitation lines, we find 
that it can fit all the high-/ CO lines observed by HIFI by sim- 
ply lowering the kinetic temperature for the inner shell from Tn 



T mb (K) 




50 100 50 100 



Fig. 7. Results of the model for IRC+10420 (red line) overlaid 
on the line profiles presented in Fig|5] as well as recent IRAM 
30m data of the /=l-0 and /=2-l transitions from Quintana- 
Lacaci et al. (in preparation) and the 12 CO J =3-2 and /=4-3 
transitions from the JCMT archive. Note that no data were found 
for the J =3-2 and /=4-3 transitions of I3 CO. For each line, 
the model was adjusted by a free corrective factor of at most 
the absolute-calibration accuracy error, which was applied to the 
given transition. For the JCMT archive data, this calibration ac- 
curacy was unclear to us so we used a conservative upper limit of 
the 25% error. For the IRAM 30m data, calibration errors of 20% 
were assumed at both 3 mm and 1 mm. The correction factor is 
indicated in the upper right corner of each box. 



= 230K to T n = 170K in order to fit the 12 CO /=16-15 line. 
This resulted in an average kinetic temperature of ~ 700 K in 
the innermost part of the CSE. We also note that the high-/ line 
intensities are fully reproducible by restricting the model to an 
outer radius of l.lxl0 17 cm, i.e. roughly the inner shell observed 
by CC07. This indicates that the outer shell around IRC+10420 
is not probed by the Herschel/HIFI CO observations. 

The result of the fitting of the HIFI 12 CO and 13 CO lines 
up to /= 16-15 from IRC+10420 is shown in Fig [7] The cor- 
responding mass-loss and kinetic temperature profiles are dis- 
played in Fig[8] We found that two different temperature dis- 
tributions can reproduce the data (Fig©. In the first case, we 
adopted the same (relatively steep) slope a = —1.2 as CC07, 
which led to a global decrease in the temperature to Tn = 170 K. 
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Fig. 8. Upper panel: Mass-loss pattern as a function of stellar 
radius derived for the model of IRC+ 10420 (Section l4~2l . Lower 
panel: Kinetic temperature profile as a function of stellar radius 
for the model of IRC+ 10420. The solid line corresponds to Tn 
= 170K and a = —1.2, while the black dashed line uses Tn = 
230 K and a = -0.8. The red dashed line illustrates the inner 
layer profile assumed in the model of CC07 (Tn = 230 K and a 
= -1.2). 



T mb /T mb (C0(6-5)) 
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Fig. 9. Histograms of the peak brightness temperatures as tabu- 
lated in Tables [3] and |H for lines detected in at least 2 sources, 
normalised by the 12 CO 7=6-5 intensity of the respective 
sources. Non-detected lines are not represented in this graph. 



However, we were also able to reproduce the observed lines we 
we considered a more gentle increase in the temperature toward 
the centre, of a = -0.8 and kept the original value of Tn - 
230 K. In both models, the highest temperature of the circum- 
stellar gas around IRC+10420 remains 5 700 K, instead of the 
values > 1000 K originally assumed by CC07 (see FiglSJ. 

The assumed 12 CO abundance was 3x10 4 (typical of O-rich 
envelopes), and the 13 CO abundance was taken from the study of 
Quintana-Lacaci et al. (2007). This translates into a 12 CO/ 13 CO 
abundance ratio of order 1/5, which is relatively high compared 
to AGBs but very close to that reported in high-mass post-red 
supergiants (see CC07, Quintana-Lacaci et al. 2007). 



5. Conclusion 

We have reported Herschel/HIFl observations of high spectral 
resolution for FIR/submm molecular lines in two red supergiant 
stars (RSGs), NML Cyg and Betelgeuse, and two yellow hyper- 
giant stars (YHGs), IRC+10420 and AFGL 2343. Tables[3]and|4] 
summarize the observational parameters derived from our data, 
which are displayed in Figs. [2]to[6]— the full spectra are shown 
in the appendix. 

As illustrated in Fig.|9j the RSGs exhibit in general more in- 
tense high-excitation lines, indicating that they harbour inner cir- 
cumstellar layers with high temperatures close to or even higher 
than 1000 K (Sections EQ and 03. In contrast, YHGs do not 
contain this central warm material, and are instead surrounded 
by detached envelopes with low densities at small distances from 
the star (Sections l3.3l and [3~4b . As such, our results show that the 
mass-loss rates in YHGs are presently weak. We also observe 
that the spectral line profiles in RSGs tend to become narrower 
with increasing excitation temperature. We show that this phe- 
nomenon indicates that the high excitation lines in RSGs orig- 
inate from gas layers at small radii that have not yet reached 
the circumstellar terminal velocity. This trend is not observed in 
YHGs, in line with the picture of a hollow envelope closer to the 
stellar photosphere. These YHGs are thus less likely than RSGs 
to contain gas still in the acceleration phase. 

In both YHGs, we detected intense emission features that 
are conspicuous in the profiles. That of IRC+10420 is relatively 
blue-shifted, at Vlsr ~ 65 km s _1 , and is more noticeable in 
certain molecules, such as H2O and NH3, although it does not 
seem to require particular excitation conditions. It was detected 
and mapped in low-/ CO Lines by CC07, who concluded that it 
comes from a condensation in the outer envelope. We also found 
a peculiar emission excess in AFGL 2343 (Section [3~4l i. at an ex- 
treme positive velocity of Vlsr ~ 125 km s _1 . There, the feature 
is clearly associated with the excitation of the lines. We specu- 
late that this shell could result from a shock interaction between 
the circumstellar envelope and nearby interestellar gas. 

We stress the particularly strong emission of water lines in 
NML Cyg and IRC+10420, which are also the only two sources 
featuring NH3 emission. On the other hand, OH is ubiquitously 
observed. It is unclear where this molecule predominantly arises 
and how it forms in the envelopes. Its emission does not seem 
to be correlated with that of H2O, although OH formation in 
evolved stars is thought to result from H2O photodissociation; 
for instance, the observed OH profiles are usually closer to those 
of mid-excitation CO lines and the OH line intensity is high in 
AFGL 2343, which shows particularly weak water lines. 

Finally, we have proposed a preliminary model to fit the 
12 CO and 13 CO lines detected in IRC+10420. The model is 
based on that developed by CC07 to explain their mm-wave in- 
terferometric mapping. We found that the original description 
of the shell around IRC+10420 by CC07 is able to reproduce 
all transitions by simply decreasing the temperature of the inner 
layers by about 30% (Fig. [8). In addition, the emission from the 
high-/ lines is found to originate only from a detached, hot shell 
close to the star, formed by stellar mass ejection at a high rate 
of ~ 3xl0~ 4 M Q yr _1 . Our calculation indicates that no hotter 
components at small radii are needed to reproduce the molecu- 
lar emission of IRC+10420, confirming that the heavy mass-loss 
ceased about 200 years ago. More detailed modelling of the CO 
lines observed in IRC+10420 and AFGL 2343 will be presented 
in a forthcoming paper, in particular to take into account the de- 
parture from spherical symmetry illustrated by high-resolution 
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maps at various wavelengths, and corroborated by the distinct Appendix A: Full HIFI spectra 

spectral features detected in the HIFI data. 
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Fig.A.l. Full Wide-Band Spectrometer spectra obtained towards NMLCyg. The setting number is indicated on the right side. 
Each spectrum is displayed on a double-side-band scale with both side-band frequency scales drawn. The spectral resolution is the 
native one (1.1 MHz), and low-order polynomial baselines have been subtracted. Lines are indicated with the respective ticks at the 
expected (velocity-corrected) frequency - lines from the USB have their ticks pointing upwards, while lines from the LSB have 
ticks pointing downwards. The spectrum in setting 16 is made of only a fraction of the data taken at this frequency, for the purpose 
of standing wave mitigation (see text for details). The obsid's numbers are given below the setting number, in the form of 1342x, 
'x' being the identifier indicated on the plot. 
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Fig. A.2. Same as Fig. lA.ll for IRC+ 10420. The spectrum in setting 16x is made of only a fraction of the data taken at this frequency, 
for the purpose of standing wave mitigation (see text for details). The tick displayed for 13 CO 7=16-15 illustrates the expected 
position of this line, despite its non-detection. 



12 



D. Teyssier et al.: HerschelfHWl observations of red supergiants and yellow hypergiants 



Tub (mK) 


1843 


1842 


1841 




: ' 1 oh a n 1/a 3/s-i/z 


ii 


.12 1 1 J 

1 CO 16-15 


500 






\ 









1 " 




1834 


1835 


1836 


Tub (mK) 


1769 


1768 





400 
200 


-200 





1760 



1761 



1762 



16 

205765 
LSB (GHz) 

USB (GHz) 

- 03b 
218422 

LSB (GHz) 



T UB (mK) 1209 



1208 



1207 



1206 



1205 




(GHz) 



1193 



1194 



1195 



1196 



1197 



LSB (GHz) 



Tub (mK) 



1165 



1164 



1163 



1162 



500 




T 



1150 



1151 



1152 



1153 



USB (GHz) 

06d 
195057 

LSB (GHz) 



T MB (mK) 1115 



1114 



1113 



1112 



1111 



Tub (mK) 




ifWilwwWi<*i'<ii%Jj*^^ 

679 680 681 

T UB (mK) 661 660 659 

150 



100 - 

50 E- 


646 



T 



CO 6-5 



647 



648 



649 



Tub (mK) 



572 



571 



570 



569 



USB (GHz) 
12f 

194541 
LSB (GHz) 

USB (GHz) 




557 



558 



559 



560 



LSB (GHz) 



Fig.A.3. Same as Fig. lA.ll for Betelgeuse. The tick displayed for NH3 illustrates the expected position of this line, despite its 
non-detection. 
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Fig.A.4. Same as Fig.|AJ]ibr AFGL2343. The ticks displayed for NH 3 , 13 CO 7=10-9, 12 CO 7=16-15 and 13 CO 7=16-15 
illustrate the expected position of these lines, despite their non-detection. 
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